Double Kelvin Wave Cascade in Superfluid Helium 
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We study the double cascade of energy and wave action in a local model of superfluid vortex fila- 
ments. The model is obtained from a truncated expansion of the 2D Local Induction Approximation 
and it is shown to support six-wave interactions. We argue that, because of the uncertainty in the 
vortex core profile, this model has the same status of validity as the traditionally used Biot-Savart 
model with cutoff, but it has advantage of being much simpler. Our minimal model leads to a wave 
kinetic equation for which we predict existence of two distinct power-law scaling in the spectrum, 
corresponding to a direct cascade of energy and an inverse one of wave action. Direct numerical 
simulations confirm the theoretical predictions in the weak turbulence regime. 



It is well known that a classical vortex filament can 
support linear waves. These were predicted by Kelvin 
more than one century ago and experimentally observed 
about 50 years ago in superfluid 4 iJe. At very low tem- 
perature, where the friction induced by normal fluid com- 
ponent can be neglected, Kelvin waves can be dissipated 
only at very high frequencies by phonon emission 
Therefore at lower frequency energy is transferred among 
different wavenumbers by the nonlinear coupling. This is 
the mechanism at the basis of the Kelvin wave cascade 
which sustains superfluid turbulence @, [I[ . 

In recent years, single vortex Kelvin wave cascade has 
attracted much theoretical [3] , numerical d, H, 0] and ex- 
perimental [|| attention. Even within the classical one- 
dimensional vortex model, different degrees of simplifi- 
cation are possible. For small amplitudes, the vortex 
configuration can be described by a two component vec- 
tor field, made of the coordinates of the vortex line in 
the plane transverse to the direction of the unperturbed 
filament. These depend on the single coordinate that 
runs along the filament. This system of equations admits 
a Hamiltonian formulation, dubbed the two-dimensional 
Biot-Savart formulation (2D-BS), see ^ below. An- 
other, more drastic, simplification is obtained by con- 
sidering local interactions only. This leads to the local 
induction approximation (LIA) which was originally de- 
rived starting from the full 3D-BS The main lim- 
itation of LIA is that it generates an integrable system 
with infinite conserved quantities, as it is equivalent to 
the nonlinear Schrodinger equation fl2| ], and therefore 
the resonant wave interactions are absent (at all orders) 
and one cannot reproduce the phenomenology of the full 
system. For this reason LIA, despite its simplicity, is of 
little help for the study of Kelvin wave turbulence. 

In this paper we consider the simplest model for vortex 
filament able to sustain a turbulent energy cascade. The 
model is obtained in the limit of small amplitudes by a 
Taylor expansion of the 2D-LIA. The truncation breaks 
the integrability of the Hamiltonian and therefore gener- 
ates a dynamical system with two inviscid invariants (en- 
ergy and wave action). For this class of systems, whose 



prototype is the two-dimensional Navier-Stokes turbu- 
lence [13} , we expect a dual cascade phenomenology in 
which one quantity flows to small scales generating a di- 
rect cascade while the other goes to larger scales produc- 
ing an inverse cascade. The possibility of a dual cascade 
scenario for Kelvin waves turbulence has been recently 
suggested [HI, [H[ but never observed. Direct numerical 
simulations at high resolution confirm the dual cascade 
picture with spectral exponents consistent with the the- 
oretical prediction based on a six-waves kinetic equation. 

At a macroscopic level, the superfluid vortex filament 
is a classical object whose dynamics is often described by 
the Biot-Savart equation (BSE) 



i< I ds x (r — s) 
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which describes the self-interaction of vortex elements. 
The quantum nature of the phenomenon is encoded in 
the discreteness of circulation k = h/m [lj|. 

The BSE dynamics of the vortex filament admits a 
Hamiltonian formulation under a simple geometrical con- 
straint: the position r of the vortex is represented in two- 
dimensional parametric form as r = (x(z), y(z), z), where 
z is a given axis. From a geometrical point of view, this 
corresponds to small perturbations with respect to the 
straight line configuration, i.e. the vortex cannot form 
folds in order to preserve the single-valuedness of the x 
and y functions. In terms of the complex canonical coor- 
dinate w{z, t) = x(z, t) + iy(z, t), BSE can be written in 
Hamiltonian form iw = SH[w]/5w* with Q 



H\w\ = — / dzidzo 
4tt J 



l + Re{w'*{z 1 )w'{z 2 )) 



\/(zi - z 2 ) 2 + \w(zi) - w(z 2 )\ 2 

where we have used the notation w'{z) — dw/dz. The 
geometrical constraint of small amplitude perturbation 
can be expressed in terms of a parameter e{z\,z 2 ) = 
\w(zi) - w(z 2 )\/\zi - z 2 \ < 1 

An enormous simplification, both for theoretical and 
numerical purposes, is obtained by means of the so called 
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local induction approximation (LIA) [11]. This approx- 
imation is justified by the observation that (TT]) is diver- 
gent as s — > r and is obtained by introducing a cutoff at 
a = |r — s| in the integral in ([1]) which represents the vor- 
tex filament radius. However, the cutoff operation does 
not take into account the distribution of the vorticity 
inside the core of the vortex filament. Indeed, let us con- 
sider linear Kelvin waves on an infinite straight filament 
with wavevector k such that h « 1. Dispersion relation 
for the frequency of such linear waves can be found ex- 
actly from the Euler equations for several special types 
of the vortex core shapes: 
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where C = —7 + In 2 + 4 for a vortex with uniform vor- 
ticity inside a cylinder of radius a [9j (where 7 is the 
Euler constant), C = —7 + In 2 for a hollow core [Icj |. 
and C = —7 — § for the BSE model {TJ) with cutoff (see 
Appendix A). From ([3]) we see that the particular shape 
of the vortex core affects the immediate next order with 
respect to the LIA term (log) and, therefore, it affects 
the leading order of the nonlinear dynamics. On the 
other hand, the choice of the vortex core shape is diffi- 
cult and uncertain, considering that the fluid description 
itself fails within the vortex core. Thus, in the follow- 
ing we will use the simplest model which arises naturally 
from the truncated expansion of the LIA Hamiltonian. 

When applied to Hamiltonian ((2|), the LIA procedure 
gives 

H[w] = 2-£- In (jj J dz^/\ + \w'{z)\ 2 = WL[w] (4) 

where I is a length of the order of the curvature radius (or 
intervortex distance when the considered vortex filament 
is a part a vortex tangle), f3 = (n/4n)]n(£/a). Here, it 
was taken into account that because a is much smaller 
than any other characteristic size in the system, (3 will be 
about the same whatever characteristic scale I we take 
in its definition. We remark that in the LIA approxima- 
tion the Hamiltonian is proportional to the vortex length 
L[w] = J dzy/l + |u>'(z)| 2 which is therefore a conserved 
quantity. The equation of motion from Q is (we set 
f3 =1/2 without loss of generality, i.e. we rescale time as 
2(3t -> t) 



w = — 



w 



(5) 



As a consequence of the invariance under phase trans- 
formations, equation ([5|) conserves also the total wave 
action (also called the kelvon number) 0] 



N\w] 



dz\w\ 



(6) 



In addition to these two conserved quantities, the 2D- 
LIA model possesses an infinite set of invariants and is 



integrable, as it is the LIA of BSE (which can be trans- 
formed into the nonlinear Schrodinger equation by the 
Hasimoto transformation, see Appendix B). Integrabil- 
ity is broken if one considers a truncated expansion of 
the Hamiltonian ([4]) in power of wave amplitude w'(z). 
Taking into account the lower order terms only one ob- 
tains: 



HexpH = H + H 1 +H 2 = 



dz ( 1 + — \w 
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Hi 



(7) 



Neglecting the constant term, the Hamiltonian can be 
written in Fourier space as 



H eX p = I uJk\wk\ z dk 



^12341^1234 SHi^wlw^W-iWi (8) 



with u> — k 2 /2, W1234 = — ^kik 2 k^k 4 and we used 
the standard notation SH(k) = 6(ki + k 2 — k 3 — k±) 
and dfci234 = dkidk 2 dk3dk 4 . For the dispersion relation 
to ~ k 2 there is no non-trivial solution to the four-wave 
resonant conditions. There is no five-wave interaction 
either, because there are no odd-degree Hamiltonians. 
The first non-vanishing process appears to be six-wave. 
Introducing a near-identity (weakly nonlinear) canonical 
transformation Wk — > Ck using the standard strategy de- 
scribed in it is possible to transform ||5J) into (See 
Appendix C) 



H r = 



2 dk 



rf&123456Ci23456^456(^) c l c 2 c 3 c 4 c 5C6 



(9) 



where C123456 = —j^kik^^k^k^kQ. This six-order inter- 
action coefficient is obtained from coupling of two fourth- 
order vertices of H 2 . It is not surprising that the resulting 
expression coincides, with the opposite sign, the interac- 
tion coefficient of H3 in ([5]). Indeed, ([S]) is an integrable 
model which implies that if we retained the next order 
too, i.e. i?3, then the resulting six- wave process would be 
nil, and the leading order would be an 8-wave process in 
this case. Thus, the existence of the six-wave process is 
a consequence of the truncation |(7j) of the Hamiltonian. 

Hamiltonian © (or equivalently ©) constitutes the 
minimal model for Kelvin wave turbulence. It possesses 
the same scaling properties as the BSE system: it con- 
serves the energy and the wave action, it gives rise to a 
six-wave system with an interaction coefficient with the 
same order of homogeneity as the one of the BSE. A 
slight further modification should be made in the time 
re-scaling factor as f3 — (k/Att), - i.e. by dropping the 
large log factor from the original definition (to have the 
correct with respect to this log in rate of the six-wave 
process) . This minimal model is much simpler than BSE 
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and it has the same degree of validity as BSE because of 
the vortex core uncertainty discussed above. 

Physical insight on Kelvin wave turbulence is obtained 
from the wave turbulence (WT) approach which yields 
the kinetic equation which describes the dynamics of the 
wave action density n k = (|cfc| 2 ). 

The dynamical equation for the variable c k can be de- 
rived from Hamiltonian ^ by the relation idck/dt = 
8H c /8c* k , and is 

i^-~U k C k = I dfc 2 3456Cfe23456C2C3C4C5Cg(S4||(fc). (10) 



Multiplying equation ([10]) by c£, subtracting the com- 
plex conjugate and averaging we arrive at 



d{c k c* k ) 
dt 



— 61m 



^23456 Cfc23456>^fc23456<5456 (&) 



(11) 

456 — \C* k C2ClC4C 5 C e ) . 

Assuming a Gaussian wave field, one can take J&23456 



where Jfc 2 3456<5| ; 23 ~ 



to the zeroth order Ji 



CO) 



which is simplified via the 
Gaussian statistics to a product of three pair correlators, 



fc23456 ' 



J, 



i°23456 = n 2 n 3 n 4 [8 k 4 (8lSl+SlSi) 

+ S k 5 (SjSl + StSi) + S h 6 (6lS 3 5 + 6 2 5 Sf) ] .(12) 



However, due to the symmetry of (7^23456 this makes the 
right hand side of the kinetic equation (fl4|) zero. To find a 
nontrivial answer we need to obtain a first order addition 



J, 



(i) 



fc23456 



to Jfc23456- To calculate J. 



(i) 

fc23456 



one takes the 



time derivative of J/c23456i using the equation of motion 
(1 1 0[) and inserting the zeroth order approximation for 
the tenth correlation function (this is similar to equation 
\12\ but a product of five pair correlators involving ten 



wavevectors) J } 



(i) 

£23456 



can then be written as 



A simple dimensional analysis of (|14[) gives 

h k ~ k 14 n 5 k (15) 

which is the same form obtained from the full BSE [4]. 

The energy flux at wavenumber k is defined as III — 

J dk'h k /ui k > which, using (fl"5"|) , becomes nj^' ~ k 11 n\. 
By requiring the existence of a range of scales in which 
the energy flux is fc-independent leads to the spectrum 



n k 



-17/5 



(16) 



A similar argument can be applied to the wave action ([5]) 



Therefore a scale 



whose flux is U k = J dk'n k i <■ 
independent flux of wave action requires a spectrum 



n k 



.-3 



(17) 



The two spectra (|16tilT[) occur in different scale ranges 
and the two cascades develops in opposite directions, as 
in the case of two-dimensional turbulence fl3| . Among 
the two conserved quantities, the largest contribution to 
energy comes from smaller scales than those that con- 
tribute to wave action (because the former contains the 
field derivatives). Therefore, according to the Fj0rtoft ar- 
gument [IH , we expect to have a direct cascade of energy 
with fc -17 / 5 spectrum at large k and an inverse cascade 
of wave action with spectrum k~ 3 at small k. 

In the following we will consider numerical simulations 
of the system ([7]) under the conditions in which a sta- 
tionary turbulent cascade develops. Energy and wave 
action are injected in the vortex filament by a white-in- 
time external forcing <fi(z,t) acting on a narrow band of 
wavenumbers around a given kf. In order to have a sta- 
tionary cascade, we need additional terms which remove 
H and N at small and large scales. The equation of 
motion obtained from (J7JI is therefore modified as 



7 (1) - Be lAuJt 

J fc23456 — JDC 



^4fc23456 



Alo 



(13) 



uj 3 - 
l 



OJ4 — 



J_ 

"3 



ojq and ^^23456 
1 1 i_ 

ri4 ns riQ 



where Alu = u k + co>2 

3 Cfc23456 n fe n 2"-3«4^5«6 

The first term of JT5|) is a fast oscillating function, its 
contribution to the integral (TTTj) decreases with z and 
is negligible at z larger than l/uj k , and as a result we 
will ignore the contribution arising from this term. The 
second term is substituted in equation pi[) .the relation 
Im(Auj) ~ — ttS(Alj) is applied because of integration 
around the pole, and the kinetic equation (fT4|) derived, 



h k = 187T / dk 



'23456 



\C k 



23456 



2 cfc23 
°456 



(A;) 



where 
n k n 2 n 3 ri4n 5 n e 



x d 456 ( w ) Jfe23456 

we have 



(14) 



introduced /fe 2 34_56 
1,1,1 1 1 i_ 

n k n 2 n 3 TI4 n 5 n 6 



1 



,/|2 



{-\yvV Ip w-aw + (j) (18) 



In (fl8|) the small scale dissipative term (with p > 1) 
physically represents the radiation of phonons (at a rate 
proportional to v) and the large scale damping term can 
be interpreted as the friction induced by normal fluid at 
a rate a. 

Assuming the spectra (|16H17[) a simple dimensional 
analysis gives the IR and UV cutoff induced by the dis- 
sipative terms. The direct cascade is removed at a scale 
k v ~ iy~ 5 /( 10 P~ 2 ) while the inverse cascade is stopped at 
k a ~ a 1 / 2 . Therefore, in an idealized realization of in- 
finite resolution one would obtain a double cascade by 
keeping kt = 0(1) and letting a — > 0. In order to have 
an extended inertial range, in finite resolution numerical 
simulations we will restrict to resolve a single cascade by 
putting kf ~ k a or kf ~ k v for the direct and inverse 
cascade respectively. 

We have developed a numerical code which integrates 
the equation of motion (|18p by means of pseudospectral 
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method for a periodic vortex filament of length 2tt with a 
resolution of M points. Linear, dissipative terms are in- 
tegrated explicitly while the nonlinear term is solved by a 
second-order Runge-Kutta time scheme. Vortex filament 
is initially a straight line (w(z, t = 0) = 0) and long time 
integration is performed until a stationary regime (indi- 
cated by the values of H and N) is reached. The ratio be- 
tween the two terms in the expansion is if i /if 2 — 20, 
confirming a posteriori the validity of the perturbative 
expansion ([7]). 

The first set of simulations is devoted to the study of 
the direct cascade. Energy fluctuations are injected at 
forcing wavenumber kf ~ 2 and the friction coefficient 
a is set in order to have k a ~ kf. Energy is removed 
at small scales by hyperviscosity of order p = 4 which 
restricts the range of dissipation on a narrow range of 
wavenumber close to k max . 



1 




1 10 100 1000 
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FIG. 1: Wave number spectrum rik for a simulation of the di- 
rect cascade in stationary conditions at resolution M = 2048. 
Forcing is restricted to a range of wavenumbers 1 < kf < 3 
and phonon dissipation is modeled with hyperviscosity of or- 
der p = 4. The line represents the kinetic equation prediction 
rife ~ fc -17 / 5 . The inset shows the spectrum compensated 
with the theoretical prediction. 

In Figure [T] we plot the wave action spectrum for the 
direct cascade run averaged over time in stationary con- 
ditions. A well developed power law spectrum very close 
to prediction p6[) is observed over more than one decade 
(see inset). This spectrum confirms the existence of non 
trivial dynamics with six-wave process for the truncated 
Hamiltonian Q. 

We remind that the direct cascade for the Hamiltonian 
([7]) was already discussed by Svistunov 0, 0] who also 
gave the dimensional prediction (|16[) . Numerical simu- 
lations with a discrete version of the Hamiltonian (J7J), 
which breaks the integrability of the original system as 
well, confirmed the validity of the dimensional prediction 
0- 

We now turn to the simulations for the inverse cascade 
regime. To obtain an inverse cascade, forcing is concen- 




FIG. 2: Wavenumber spectrum nt for a simulation of the in- 
verse cascade in stationary conditions at resolution M = 2048. 
Forcing is restricted to a range of wavenumbers around kf = 
300 and phonon dissipation is modeled with hyperviscosity of 
order p = 4. The line represents the kinetic equation predic- 
tion ~ k~ 3 . The inset shows the spectrum compensated 
with the theoretical prediction. 



trated at small scales, here kf = 683. In order to avoid 
finite size effects and accumulation at the largest scale 
[l9j |. the friction coefficient is chosen in such a way that 
wave action is removed at a scale k a ~ 10. Figure [5] 
shows the spectrum for this inverse cascade in stationary 
conditions. In the compensated plot, a small deviation 
from the power-law scaling at small scales, probably due 
to the presence of forcing, is observed. Nevertheless, a 
clear scaling compatible with the dimensional analysis of 
the kinetic equation is observed over about a decade. 

In summary, we have introduced a minimal model for 
Kelvin wave turbulence, and presented an argument why 
this model should be preferential over BSE. Namely, this 
model is much simpler than BSE and it has the same 
degree of validity due to the vortex core shape sensitivity 
and uncertainty. We have used our model for numerical 
simulations of the direct and the inverse cascades and 
found spectra which are in very good agreement with the 
predictions of the WT theory. 

I. APPENDIX A - INTERACTION 
COEFFICIENTS OF 2D-BS 

In this Appendix we extend the work of Kozik and 
Svistunov on the Kelvin wave cascade (KS04) 0]. They 
consider the full Biot-Savart Hamiltonian © in 2D, and 
simplify the denominator by Taylor expansion. The cri- 
terion for Kelvin- wave turbulence is that wave amplitude 
is small compared to wavelength, this is formulated as: 

\w{z x ) -w(z 2 )\ 
e(zi,z 2 ) = j j < 1. (19) 

\Zl — Zi\ 
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KS04 find the Biot-Savart Hamiltonian ([2]) expanded 
in powers e (H = H + Hi + H 2 + H 3 , here H is just a 
number and is ignored) is represented as: 



Hi = 
H 2 = 

H, = 



k f dzidz-i 



■$tt J \z\ - z 2 
/.' f dz\dz2 



2Re I w *{zi)w (z 2 ) 



32-71 J \zi — z 2 \ 
dz\dz 2 



3e 4 - Ae 2 Re 



647T J \zi - z 2 \ 



6e 4 Re 



(w *(zi)w (z 2 ) 
(w *{zi)w (z 2 ) 



(20) 
5e 6 l 



By taking the Fourier transformation w(z) = 
J dkw k e lkz of (|20|) , invoking a cutoff at \zi — z 2 \ = a be- 
cause of the singularity present in the Biot-Savart Hamil- 
tonian @ as \zi — z 2 \ ~ > 0, KS04 derived the coefficients 
of H\ , iJ2 and i/3 in terms of cosines in Fourier space [|[ . 
Kozik and Svistunov however, left their work unfinished. 
We extend the work of KS04 by integrating the Fourier 
representations of Hi, H 2 and H 3 using integration by 
parts, and applying the following cosine identity [2(]| 



cos(t) 



dt 



-7 - ln(a) - / — ^-dt 



-7 - ln(a) - ^ ' 



tri 2fc(2fc)! 
= -7 - ln(a) + 0{a 2 ) (21) 

Neglecting terms of order ~ a 2 and higher we calculate 
the frequency and interaction coefficients of equations 

mm 



Hi 
H 2 
H 3 
where 



Lu k \w k \ 2 dk, (22) 
dki 23 4Wi 2 34Wiw 2 wlwlSH(k), (23) 

^123456^123456^1 ^2 IU5 lUg (fc) (24) 





3" 




~ 7 ~2 





ir, 



1234 



Ci 



—k 2 

Air 



——kik 2 k 3 k 4 

647T 

-F1234, 

K 



1 + 47 - 4 In 



1 



23456 



1287T 
Gi23456- 



kik 2 k 3 k 4 k 5 k 6 



k e ga 
1 - 47 + 4 In 



(25) 



(26) 



kr.a a 



(27) 



We use the notation that k fr is the mean value of 
wavenumbers, 7 = 0.5772 ... is the Euler Constant and 
F1234 and G123456 are logarithmic terms of order one 



shown below 

-F1234 



_/V 4 

167r [ 6 E^I ln 

Nei 



N 

fceff 



N 



Ekjki In , 

NeJ 



(28) 



G 123 456 = 7^[ 3 E^2^1n(-^ 



^ 6 , ( N 
5 2. 720^ — 



I = {-li], 



r4i r3i r43i r4 
I J- [2J- 



(29) 
(30) 



: { {H, [l]>-[ 43 ]>-[2]} l= 4 !i=1 > 

{[ 3 ],[l] ; -[ 43 ],-i} i=3 , j=1 . 
{[ 3 ]j2],-[ 43 ],-[?]} 



i=3,j=2 



K 



{[%[2] 



[ 4 ]} 



[IV 



i=4j=2 J 



(31) 



r56i 



EL 



in.-ri.H.B 8 ]. 



(32) 



r65i 



[iU 4 L-[%[ 46 Ui],-[ 46 ],- 



r456i 



[l 6 ],[23] 



L3J, L 3 J, L13J 



[ 5 ],[ 45 ],[?],-[f], 



[13]; [3]) [3 s ]' 



(33) 

The notation used for the logarithmic terms is: for N £ 
K = ±[^ 9 ] then the corresponding term in G123456 is 

— log 



3k (k p 
16tt b 



24 



As you can see, the logarithmic terms F1234 and G123456 
are of order ~ 1 and are extremely messy and complex. 
Applying the parametrization of Zakharov and Schulman 
[211 ] does not help in this case (see Appendix C). We 
know that LIA is fully integrable, i.e. that after apply- 
ing the canonical transformation all terms of order ln(fca) 
will cancel. However, the presence of the order one ad- 
ditions to the dispersion and interaction coefficients, will 
break this integrability and result in some six- wave inter- 
actions. Due to us being unable to simplify the expres- 
sion after the canonical transformation because of the 
non-canceling order one terms, we make the assumption 
of the interaction coefficients being products of wavenum- 
bers, this is equivalent to the truncated LIA model sug- 
gested in this letter. 



II. APPENDIX B - INTEGRABILITY OF 2D-LIA 

We show the details of the derivation that the LIA of 
BSE yields the 2D-LIA model ©, i.e. that the cutoff 
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operation commutes with making the 2D reduction. In 
the view of integrability of the original LIA, this amounts 
in a proof of integrability of the 2D-LIA model ([5|) . 
The LIA representation of BSE can be written as 



/3r' x r" 



(34) 



where the notation for the differentiation ' is d/dl 

— 1/2 

where / = (l + |w| 2 ) is the arc length. The two- 
dimensional representation of a vortex line can be rep- 
resented by a vector r = zz + w. The vector w = 
(x(z), y(z)) being a function of z, orientated in the xy- 
plane. 

Applying the chain rule to rewrite all derivatives to be 
with respect to z (from this point on ' will refer to d/dz), 
LIA can be represented as 



r = /3 (l + |w'| 2 ) 3/2 (z x w" + w' x w") 



(35) 



With a little geometrical intuition, one can show w = 
r — (r • z ) (z + w' ) [2| . Both w' and w" are perpendicular 
to z direction, thus, one can represent w' x w" = ((w' x 
w") • z)z = Az. Then equation (|35|) reduces to 

w = p (1 + |w'| 2 )~ 3/2 [z x w" - Aw'] . (36) 

We will now show that equation ([5|) is equivalent to 
equation (f36|) . First we must change our representation 
of w from a complex variable w(z) = x(z)+iy(z) to vector 
notation w = [x{z),y{z)). Equation §5§ is equivalent to 



1„ d / w' 
w = —z X — — I = 
2 dz I + 



(37) 



Expanding, keeping track of (3 and applying the vector 
identity (a x b) x c = (c • a)b — (c • b)a, equation (f37|) 
can be re-written as 



w = /3(l + |w' 



-3/2 



[z x w" - Aw'] 



(38) 



This is exactly the same result as Equation (J36J). The 
2D-LIA model © is equivalent to the LIA of BSE (J34J) 
and so the 2D-LIA model is indeed integrable. 



III. APPENDIX C - CANONICAL 
TRANSFORMATION 

In Wave Turbulence, the near-identity transformation 
allows one to eliminate "unnecessary" lower orders of 
nonlinearity in the system if corresponding order of the 
wave interaction process is nil [17| . For example, if there 
is no three-wave resonances, then one can eliminate the 
cubic Hamiltonian. (The quadratic Hamiltonian corre- 
sponding to the linear dynamics, of course, stays). This 
process can be repeated recursively, in a way similar to 
the KAM theory, until the lowest order of the non-trivial 



resonances is reached. If no such resonances appear in 
any order, one has an integrable system. 

In our case, there is no four-wave resonances (there are 
no non-trivial solution for the resonance conditions in for 
u> ~ k x if x > 1). There is also no five- wave resonances 
because the original terms in the Hamiltonian are of the 
even orders. However, there are nontrivial solutions of 
the six-wave resonant conditions. Thus, one can use the 
near-identity transformation to convert our system into 
the one with the lowest order interaction Hamiltonian to 
be of degree six. Let us do this. 

A trick for finding a shortcut derivation of such a 
transformation was found in [ijj • It relies on the fact 
that the time evolution operator is a canonical transfor- 
mation. Taking the Taylor expansion of w(k, t) around 
w(k, 0) = c(fc,0) we get a desired transformation, that 
is by its derivation, canonical. The coefficients of each 
term can be calculated from an auxiliary Hamiltonian 
H aU x- Similar procedure was done in Appendix A3 of 
[17] ] to eliminate the cubic Hamiltonian in cases when 
the three-wave interaction is nil, and here we apply a 
similar approach to eliminate the quadric Hamiltonian. 
The transformation is represented as 



Wfc = c (M) + t( 9c(M) 



Of 



Z 2 
2 



d 2 c(k,t) 
dt 2 



(39) 

The transformation is canonical for all t, so for simplicity 
we set t = 1. The coefficients of (|39|) can be calculated 
from the followinf formulae, 



f de(Je,t) 
\ dt 
d 2 c(k,t) 
dt 2 



.SH„ 



t=o 



t=o 



5c* 
. d 6H av 
' dt 5c* 



(40) 



Due to the original Hamiltonian ([7]) H exp having U(l) 
symmetry we have no odd wave-interactions, this sim- 
plifies the canonical transformation (|3T)|) greatly, because 
the absence of odd interaction in H exp automatically fixes 
the arbitrary odd interaction coefficients within the aux- 
iliary Hamiltonian H aux . Transformation (|39p reduces 
to 



W k = C k - - J dk 2 34Wk234C2C3C45^(k) 

-3i J ^23456^23456^456 ( fc ) c 2 c 3 c 4 c 5C6 

+ 1 J dk 234567 (w k 743Wt 623 5f 5 (k)5 7 2 l(k) 

-2WWVF7356 X 5%(k)5ll(k))c* 2 C* 3 C4C 5 C 6 .{41) 



Here we have used tildes to represent interaction coeffi- 
cients of the auxiliary Hamiltonian H aux . To eliminate 
the nonresonant interactions we substitute transforma- 
tion (|4"Tj) into Hamiltonian ©, this will yield a new rep- 
resentation in variable Ck for Hamiltonian (151) where the 



7 



nonresonant terms (more specifically the four-wave inter- 
action terms) will involve both W1234 and W1234. Arbi- 
trariness of W1234 implies that we can select this to elim- 
inate the total four-wave nonresonant interaction term, 
H 2 present within the Hamiltonian. In our case, 



W?234 = 



4iW 



1234 



U>1 + UJ 2 — UJ3 — LO4 



(42) 



This choice is valid as the denominator will not van- 
ish due to the nonresonance of four-wave interactions. 
Hamiltonian ([8]) expresses in variable c k , H c becomes 



H c = J u k c k c* k dk 

+ / ( C123456 - « (Wi + W 2 + W 3 - W 4 - W 5 - UJ 6 ) 



X C123456 <$456 ( k )c*iC* 2 C* 3 C 4 C 5 C 6 dk 123 4 56 , 



(43) 



C123456 is the arbitrary six-wave interaction coefficient 
from the auxiliary Hamiltonian. This term does not con- 
tribute to the six-wave resonant dynamics as the factor 
in front will vanish on the resonant manifold that ap- 
pears in the kinetic equation. However, we can select 
i (uji + u 2 + lo 3 - lo 4 - uj b - uj 6 ) C123456 to equal C123456 
off the resonant manifold. This enables us to write H c 
as equation ([5]), where the explicit formula for C123456 
stemming from the transformation is 



C 



123456 



6 

E 

i,j,k—l p,q,r— 4 
i^j^k^i p^q^r^p 



E 



+ 



ppij ^y q-\-r — kkqr 



Uk—Uq- 



(44) 



Zakharov and Schulman discovered a parametrisation 
PH for the six- wave resonant condition with ui k ~ k 2 , 
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(45) 



This parametrisation allows us to explicitly calculate 
C123456 on the resonant manifold. This is important 
beecause the wave kinetics take place on this manifold, 
that corresponds to the delta functions within the ki- 
netic equation. When this parametrisation is used with 
equation (|44|) and W1234 = —^kik 2 k 3 k4 we find that 
the resonant six-wave interaction coefficient simplifies to 
C123456 = ^\\kik 2 k 3 k4k 5 kQ. Note, that this is indeed 
the identical to the next term, H3 in the LIA expansion 
with opposite sign. 
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